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. It is well known that the dynamics of three point vortices moving in an ideal fluid in the plane can be 
expressed in Hamiltonian form, where the resulting equations of motion are completely integrable in the 
sense of Liouville and Arnold. The focus of this investigation is on the persistence of regular behavior 
(especially periodic motion) associated to completely integrable systems for certain (admissible) kinds 
of Hamiltonian perturbations of the three vortex system in a plane. After a brief survey of the dynamics 
of the integrable planar three vortex system, it is shown that the admissible class of perturbed systems 
is broad enough to include three vortices in a half-plane, three coaxial slender vortex rings in three- 
space, and 'restricted' four vortex dynamics in a plane. Included are two basic categories of results 
for admissible perturbations: (i) general theorems for the persistence of invariant tori and periodic 
orbits using Kolmogorov-Arnold-Moser and Poincare-Birkhoff type arguments; and (ii) more specific 
and quantitative conclusions of a classical perturbation theory nature guaranteeing the existence of 
periodic orbits of the perturbed system close to cycles of the unperturbed system, which occur in 
abundance near centers. In addition, several numerical simulations are provided to illustrate the 
validity of the theorems as well as indicating their limitations as manifested by transitions to chaotic 
dynamics. 

Physics and Astronomy Classification Scheme: 47.10.Df, 47.10.Fg, 47.15.ki, 47.20.Ky, 47.32. C-, 
47.32.cb, 47.52. +j 

I. INTRODUCTION 

Although virtually all research on vortex dominated fluid flows has its roots in the seminal work of 
Helmholtz |29| (c/. |38|). specific advances in the dynamics of point vortices moving in an ideal (= inviscid, 
incompressible) fluid in a planar region , which wc shall refer to as n vortex dynamics or the n vortex problem, 
can be traced back to the pioneering work of Kirchhoff 33 and Grobli j^Hl ■ Kirchhoff was the first to describe 
the Hamiltonian structure of the n vortex problem, which he used to derive some fundamental integrals of 
the motion, while Grobli conducted a detailed analysis of the three vortex problem that included what 
was essentially a proof of the integrability by quadratures of the problem without employing Hamiltonian 
formalism, although he did not give a complete description of the dynamics. About twenty years after 
Grobli's remarkable work, Goryachev |27| took up the problem again, and was able to obtain new insights 
concerning the dynamics of point vortices. Shortly thereafter, Poincare |47| put his own imprimatur on 
vortex dynamics, just as he did in so many other fields of research. 

Some seventy years after the pioneering work of Kirchhoff and Grobli, Synge |52| was able - using 
trilinear coordinates - to fill many of the gaps in the dynamical picture of the three vortex problem left 
by earlier studies. Among Synge's most important contributions were the derivation of integrals and the 
characterization of the critical points in terms of trilinear coordinates, and the identification of a single 
parameter - involving the sum of the product of the vortex strengths - that distinguishes three distinct types 
of qualitative dynamics (the elliptic, hyperbolic, and parabolic types) depending on whether this parameter 
is positive, negative or zero. Twenty years later, Aref \I\ E] used the Hamiltonian based investigations of 
Novikov to rediscover the advantages of studying the three vortex problem in the context of trilinear 
coordinates. In the process, Aref identified additional properties of three vortex dynamics, and initiated 
research on chaos in point vortex dynamics and its relation to turbulent flows. 

Subsequently, Tavantzis & Ting |53j . taking Synge's approach as their point of departure, derived a 
new constant of motion in trilinear coordinates, which they used together with some classical perturbation 



techniques to nearly complete the characterization of three point vortex dynamics. Using linear analysis, they 
determined the stability of all isolated stationary points in the trilinear plane, and showed that expanding and 
contracting configurations are, respectively, stable and unstable, which implies that contracting similarity 
solutions and the eventual collision of three vortices is unstable. More recently, Ting et q/.|56j employed 
techniques from nonlinear stability analysis to supply the few missing details in the dynamics; in particular, 
they showed that orbits starting just off the contracting configuration branch of the singular curve of critical 
points in the parabolic case are ultimately attracted to the expanding branch of this curve. 

The Hamiltonian approach introduced by Kirchhoff |33], further developed by Lin and perfected for 
point vortex problems by Novikov |45| has proven to be very useful in vortex dynamics research. Although 
not essential for solving the three vortex problem - as demonstrated in [2H|, E2], EH], and \5Q - presumably 
one could use the integrals in involution to reduce it to a solvable one-degree-of-freedom Hamiltonian system 
along the lines indicated by Borisov and his collaborators in such papers as E]> CB]> and E3- ^ n ^ ne other 
hand, the symplectic structure underlying Hamiltonian dynamics has proven to be extraordinarily effective 
in resolving a wide variety of vortex problems such as the formulation of point vortex dynamics on the sphere 
by Bogomolov , the proof of complete integrability of the three vortex problem on the sphere by Kidambi 
& Newton |M] (see also [T^], \T7\, [35] and 03]), verification of non-integrability of the general n vortex 
problem and (n — 1) coaxial vortex ring problem for n > 3 by Bagrets & Bagrets jS] (see also Ziglin |58|h 
and several other important results such as in |3], 0, ES], [23, EH], H§|, and |57| . 

It is in the Hamiltonian perturbation of integrablc point vortex dynamics where the Hamiltonian approach 
has proven to be particularly useful. This is largely due to the availability of two of the most important 
results in finite-dimensional symplectic dynamics: Kolmogorov-Arnold-Moser (KAM) theory (such as in 0], 
[3T], and 021) and Poincare-Birkhoff (PB) theory and its extensions and variants (see e.g. [7], |1U| . 
EU, EH, ED> E2, EH, E3, E3> EDI, EB, E3, ESI, and [H]). Examples of such applications manifold: 
Khanin [22 used a KAM theory inspired method to show that there exist subsets of initial configurations 
for the (non-integrable) four vortex problem in the plane leading to regular (integrable like) motion; namely, 
quasiperiodic orbits on (deformed) invariant (KAM) tori (cf. Celletti & Falcolini [21). A KAM theory 
based argument combined with a deft application of Jacobi canonical transformations enabled Lim |.'i'J| to 
prove the existence of quasiperiodic flow regimes in lattice vortex systems. Blackmore & Knio [jj], using 
the fact that slender coaxial vortex ring dynamics can be viewed as a perturbation of point vortex motion, 
employed an innovative KAM theory type result to prove the persistence of KAM tori and periodic orbits for 
an ample set of initial positions of three rings sufficiently close to one another when the rings having vortex 
strengths of the same sign. Later, Blackmore et al. employing an analogous KAM theory approach in 
concert with a novel extension of PB theory, generalized this to any finite number of coaxial vortex rings 
with strengths of the same sign. Similar results were obtained by Blackmore & Champanerkar using the 
same type of approach for any finite number of point vortices in a plane or half-plane. In addition, Ting & 
Blackmore [HE] and Blackmore et al. ^2] employed analogous ideas, together with trilinear coordinates, to 
prove the persistence of regular flow regimes for three coaxial vortex rings and three vortices in a half-plane, 
respectively, when the vortex strengths are of differing signs. 

In this paper we shall, after a brief description of three vortex dynamics, assemble and extend some 
of our recent results on the persistence of quasiperiodic flows on KAM tori and periodic orbits for certain 
types of non-integrable Hamiltonian perturbations of three vortex dynamics. Although we shall include 
some interesting - fundamentally qualitative - results on the existence of quasiperiodic and periodic orbits 
for perturbations of three vortex dynamics, we intend to focus our attention on more quantitative behavior 
concerning the persistence of periodic orbits that are, in an appropriate sense, near periodic orbits of the 
unperturbed completely integrable system. We shall obtain our results using a judicious combination of 
Hamiltonian techniques and more classical perturbation methods that are closely linked with the represen- 
tation of the unperturbed dynamics in trilinear coordinates. 

The equations of motion for three vortex dynamics in both Hamiltonian form and planar trilinear coor- 
dinates are presented in Section II, along with a brief review of the trilinear phase portrait that emphasizes 
those dynamical features that exploited extensively in the sequel. Next, in Section III, we describe specific 
types of perturbations that are subsumed by the Hamiltonian perturbations for which our main conclusions 
are proven in subsequent sections. The specific perturbations considered are three vortices in a half-plane, a 
'restricted' four vortex problem in the plane, and three coaxial vortex rings. We prove qualitative theorems 
in Section IV on the persistence of quasiperiodic flows on KAM tori interspersed with periodic orbits for 
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Hamiltonian perturbations of three vortex dynamics subject to hypotheses satisfied by the examples in the 
preceding section. A more classical perturbation approach is used in Section V to show that the pertur- 
bations of type introduced in Section IV have the property of having periodic orbits that are close - in an 
appropriate coordinate system - to certain cycles of the three vortex system. 

In the penultimate part of the paper, Section VI, we present a variety of numerical simulations that 
illustrate the persistence of regularity associated with the three vortex problem under perturbations satisfying 
the properties in our main results. The properties in question provide sufficient conditions on the initial 
configurations of point vortices or slender coaxial vortex rings for the existence of invariant tori and periodic 
orbits, and for periodic orbits of the perturbed system that are close to those of the unperturbed system. A 
number of simulations are also included to show how regularity breaks down - signalled by the emergence 
of chaos - as the limitations imposed by the hypotheses in our theorems are exceeded. Finally, in Section 
VII, we summarize and distill the most important conclusions in the paper, and indicate some promising 
directions for related future research. 

II. UNPERTURBED DYNAMICS 

In this section we describe the equations of motion of point vortices in an ideal fluid in the plane, and 
give a rather complete characterization of the dynamics for the integrable three vortex system. Knowing 
the dynamics in the three vortex case shall prove quite useful in our subsequent analysis and description of 
periodic orbits for perturbed systems that are perturbations of cycles of the unperturbed system. 

A. Governing equations 

We begin with the general problem of n point vortices of respective nonzero strengths Ti, T n moving 
in an ideal fluid in the (complex) plane C (= K 2 ) and located at the positions z k — x k + iy k , 1 < k < n, 
respectively. The equations of motion in complex form are 

n 

h = i ^2 k j ( z j ~ z ^y x > ( i - k - n ) ( i ) 

where the overbar indicates the complex conjugate, the overdot denotes differentiation with respect to t, and 
Kj := Tj/2-K, 1 < j < n, which can be recast as the complex Hamiltonian equation 

K*i = 2id z H , (2) 

where k := (k\, K n ), z := (zi,...,z n ), * denotes the usual (Cartesian product ring) product in C", 
d z H := (d Zl H a , ...,d Zn II ), and the Hamiltonian function is given as 

H := - ^2 KjKfc log \zj - z k \ . (3) 

l<j<k<n 

We note here that to guarantee smoothness of the system (1), the phase space needs to be defined as 

C% := {z e C" : Zj ? z k Vj + k} , 
and we introduce the following notation that will prove useful in the sequel: 

n 

:= Kj, &^ := KjKk- 

3=1 \<j<k<n 

The governing equations can be converted into the more familiar real Hamiltonian form in M. 2n (or more 
properly R 2 ^) 

x k = K^dy^v = {H ,x k } , y k = -K~ k 1 d Xk H = {H ,y k } , (1 < k < n) (4) 

where is in the usual representation in real coordinates, and the (nonstandard) Poisson bracket is 
defined as 

Sf 1=Y" -i ( d J_ dff df dg 

Vd2/ fe fe fe dx k dy k 
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It is easy to see that 

n n 
k=l k=l 

representing the total energy, the linear momentum (impulse) and the angular momentum, respectively, are 
integrals of the system (1) (— (2) = (4)). From these we can construct the following three functionally 
independent, real constants of motion that are in involution: 

H , K, L := (<HJ) 2 + (JJ) 2 , (6) 

where £H and ^denote the real and imaginary part, respectively, of a complex number, and it is easy to verify 
the involutivity of these integrals by computing that 

{H o ,K} = {Ho,L} = {K,L} = 0. (7) 

Accordingly the dynamics of three point vortices in a plane is completely integrable in the sense of Liouville 
and Arnold, which we shall denote as LA-integrable. Whence it follows from Liouville- Arnold theory and the 
Poincare-Birkhoff fixed point theorem and its extensions (see e.g. [5]. |7] P^ l - P^ l - P^l- P H - W \ - P^ l - P^ l - 
|26) . |3U| . |31|. |44|. and 021) that three vortex dynamics is quasiperiodic, confined to invariant 3-tori and 
exhibits periodic orbits of arbitrarily large periods for all combinations of nonvanishing vortex strengths. 

B. Dynamics in trilinear coordinates 

Three vortex dynamics, which we summarize here, can be described in a very efficient manner using 
trilinear coordinates, which have their roots in algebraic geometry. Our approach follows that of Synge |52| . 
Tavantzis & Ting [SH], and Ting et al. [HE]) and we refer the reader to these sources and Aref 0Ej for further 
details. We begin with a dynamic formulation introduced by Grobli [2B] in terms of the sides of the (possibly 
degenerate) triangular configuration denoted as R\ := \z% — Za\, R2 '■= \z\ — 23 1, and R% := \z\ — z%\, which 
can be expressed as 

R1R1 = 2Ak! (R2 1 - R3 1 ) , 

R 2 R 2 = 2Ak 2 (R- 1 - R- 1 ) , (8) 
#3^3 = 2Ak 3 (R^ 1 - R^ 1 ) , 

where A denotes the oriented area of the configuration defined to be positive (negative) for a counterclock- 
wise (clockwise) arrangement of the vertices ordered as z\,Zi,Zz- Naturally the points R = [R\, R2, R3) 
in Euclidean three-space that comprise the domain of (8) must be confined to the first octant (and also 
avoid or compensate for singularities on the bounding coordinate planes). It is interesting to note that the 
stationary (fixed) points of (8) correspond to equilateral configurations, or collinear configurations satisfying 
the additional condition A = 0, where the extra condition in the collinear case is necessitated by the singular 
nature of the vector field at such points. Observe that the systems (I), (4) and (8) are invariant under the 
transformation t — > —t, 

3 3 

(I < j < 3), hence we may assume without loss of generality that 

Ki > K2 > and k 2 > K3, (9) 

which we shall do throughout the sequel. 

The triangle inequality imposes additional restrictions on the domain of (8). We shall define this restricted 
domain in accordance with the approach of Synge |52j (and Tavantzis & Ting |53j). wherein one can find 
some excellent figures representing the concepts described below, by first introducing the equilateral triangle 
T := P1P2P3 of height one determined by intersecting the first octant with the plane V defined as the locus 
of R 1 +R 2 +R 3 = v/2/3, where Pi := L/2/3, 0, 0) , P 2 := (o, 0) , and P 3 := (o, 0, y/2/tj . Observe 

that the first octant can be viewed as the fiber bundle over T with the rays emanating from the origin (but 
not including the origin) as fibers. It is easy to verify from the triangle inequality that the proper base 
space of the bundle defining the admissible points V for (8) is the subtriangle T = Q1Q2Q3, where Qi is the 
midpoint of the edge P2P3, Q2 is the midpoint of the edge P1P3, and Q 3 is the midpoint of the edge P1P2, 
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so that Qi := ((3, yJTJ&, y/TJ^j , Q 2 := (v^VM, VX^) . and Q3 : = (\/V6, a/Vm). Wc note that the 
bundle projection p : I? — > T is just the radial projection on T. 

Points in the triangle T can be assigned trilinear coordinates x = (#1, £2, £3) defined as 



xj := i?^ (i?i + R 2 + R 3 ) 



(1 < 3 < 3) 



(10) 



which we note are not independent inasmuch as they satisfy x\ + x 2 + x 3 = 1, and represent distance from 
the edges of T ; in particular. x\, x 2 . and x 3 are, respectively, the distances of x from the edge's P 2 P 3 , P1P3, 
and P\P 2 - The points in the planar triangle T can also be conveniently represented in terms of the following 
coordinates: 

a := (l/Vt) (x 2 - Xi) , /? 1=0:3, (11) 

with 'inverse' transformation 



Xl = (1/2) (l - /3 - aVi) , x 2 = (1/2) (l-0 + aVty , x 3 = (3. 
The system (8) projected on T assumes the following form in trilinear coordinates 

±1 = A {kiXi (xj - xj) - xi [kiXi {x% - xl) + k 2 x 2 (xj - xj) + k 3 x 3 (xl - x\)] } , 
x 2 = A {k 2 x 2 (xl - xj) - x 2 [kiXi (xj - xl) + k 2 x 2 (xj - xl) + n 3 x 3 (xl - xj)] } , 



(12) 



x 3 = A {k 3 x 3 (xl - xl) - x 3 [mxi (xl - xl) + k 2 x 2 (x\ - xf) + k 3 x 3 (a 



xl) 



where 



A := 2A 



(RiR 2 R 3 ) 1 (Pi + R 2 + R 3 ) 



(13) 



(14) 



The dimensionality of (12) can be reduced by one by recasting the system in a, /3 - coordinates as 

a = A { (l/8V3) [/S2(l - 3/3 - aV3) (l - {fi - aV3) 2 ) + 

m{\ - 3/3 + aV3) (l - (/? + oV3) 2 ) - asj , 
/? = A { K3 V3a/3 (1 - /?) - (/3/8) s} , 



(15) 



where 



S:= j-Ki (l-3/3 + aV3) l-(/3 + a\/3) 



k 2 (l - 3/3 - aV3) 1 - (/3 - aV3) 



+ K38V3. 



a/3 (1-/3) J 



(16) 



and the factor A can be absorbed into a rescaling of time in order to simplify the description of the phase 
portrait on T. 

The system (13) or (15) actually includes two possible orientations of the triangular configurations of 
three vortices, so along with the phase portrait on T, we must include another copy of T, which we denote as 
T* , having the opposite orientation. More precisely, the integral curves on T* are precisely those of T except 
they have the opposite orientation with respect to t. By gluing T and T* (and their vector fields) smoothly 
together along their corresponding sides to obtain the double M := D(T) of T, which is a 2-sphere, we obtain 
a vector field X on M by joining the oriented field corresponding to (13) or (15). Note that at this point, 
we can reformulate the bundle structure described above for V in the 'doubled' form p : V — > M = D(T). 
As for fixed points of the flow, we always have, in terms of trilinear coordinates, E := (1/3, 1/3, 1/3) on 
T and its copy E* on the T* half of M corresponding to an equilateral configuration, and the common 
vertices Qi = Q\, Q 2 = Q 2 and Q 3 — QI corresponding to points where a pair of vortices coincide, as well 
as additional points that we shall describe in what follows. 

A local stability analysis at E (or E*) shows that there are three distinct types of dynamical behavior 
corresponding to the value of := .ft^: These are the elliptic case when > 0, the hyperbolic case 
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when < 0, and the exceptional parabolic case when = 0. When the system is elliptic, E (and E*) 
is a center; E (and E*) is a saddle point in the hyperbolic case; and when the system is parabolic, E (and 
E*) is no longer isolated - it lies on the curve of fixed points C on T (with copy C* on T*) defined as 

One of the most remarkable and useful features of the trilinear coordinate based projection of the system 
(8) on M is what amounts to essentially a unique path lifting property (cf. for integral curves on M, 

except along C for parabolic systems. For future reference, we shall refer to this property as the unique 
integral path lifting principle for the bundle p : D — ► M = D(T). More precisely, the initial points in T> 
for integral curves of (8) establishes a bijective correspondence via the projection p with the integral curves 
of (12) or (14) on M, except along C in the parabolic case. This bijective correspondence can be readily 
established using the integrals of motion of the systems. As for these flow invariants, employing the constants 
of motion Hq and K for (4), it is easy to verify that 

«i X logi?i + k^ 1 logi? 2 + logi? 3 and n^ 1 Rj + k^ x R\ + k^ x R\ (17) 

are first integrals for (8). These integrals were cleverly combined by Tavantzis & Ting 53 to obtain the 
following invariants for (13): 

( 3 \ / 3 \ ( K i K 2'W- s<2) ) 

XX 1 **) (lK /Kfc J . ( 18 ) 

when £( 2 ) ^ 0, and 

I := (zi/*3) 2K2 (x 2 /x 3 ) 2Kl (19) 

for the parabolic case. The above integrals can be combined (see ^2] an d to produce a an invariant 
valid for all cases, namely 

We note that we shall find the invariant / to be particularly useful for our numerical simulations of pertur- 
bations of three vortex dynamics in the sequel. In particular, the extent to which I differs from a constant 
value for a perturbed system is a good indication of the extent to which the dynamics of the systems diverges 
from the LA-integrable motion of the three vortex problem. Any of the above constants of motion can be 
used to prove the almost unique path lifting property (cf. Synge |52j). 

Before embarking on a condensed - but reasonably complete - trilinear coordinate based description of 
the dynamics in the elliptic, hyperbolic and parabolic cases, certain aspects of the compelling effectiveness 
of the trilinear approach almost demand comment. We first observe that (8) - and so also (13) - incorporates 
the Lie group of symmetries of the original system (1) - which is the planar Euclidean group. As one would 
expect to be able to use the symmetries in a standard symplectic fashion via the associated momentum map 
to obtain a Hamiltonian restriction of the original system (4) having just one degree of freedom for the three 
vortex problem, there is the obvious suggestion of a strong link between the symplectic reduction and the 
system (13). It appears that such a link, which shall have no need of in the sequel, has yet to be definitively 
established. 

Elliptic case (8.^ > 0) 

As indicated above, in the typical elliptic case, which is illustrated for T in Fig. 1, both the stationary 
point E and its (opposite orientation) copy E* are centers with index +1. It can be easily shown that the 
fixed points Q±, Q2 and Q3 are also centers. In addition, there are three other fixed points of (13) - one on 
the interior of each side of T - denoted as Q4, Q$ and Qe in Fig. 1, and each of these additional stationary 
points can readily be shown to be a saddle points of index —1. Notice that these observations are consistent 
with the Poincare-Hopf Index theorem (see e.g. |31|) , inasmuch as a simple calculation yields 

5(+l) + 3(-l) = 2 = X (M) = X (D(T)) , 
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where x{M) is the Euler-Poincare characteristic of the sphere M, which is two. In this case, we see that 
the global phase portrait of (15) is completely determined, and this leads to an exhaustive determination 
of three vortex dynamics in the elliptic case owing to the unique integral path lifting principle mentioned 
above. 




°3 

Figure 1: Elliptic type with k\ = 2, K2 = 1, and K3 = 1/2. 

Hyperbolic case (A^< 0) 

An example of the dynamics for the hyperbolic case is shown in Fig. 2. For all hyperbolic cases, the 
points E and E* associated to equilateral configurations are stationary saddle points of index —1. When it 
comes to the fixed points at the vertices of the triangle T (glued to those of T*), the point Q3 is always a 
center of index +1 for all possible admissible vortex strengths. Any other fixed points, which of course must 
include the other vertices of T as well other points on the open edges of the triangle can have a variety of 
natures depending on certain algebraic criteria (see For example, in the case shown in Fig. 2, Qi and 

Q2 are both centers, while the stationary points Q4, Q5 and Qe are, respectively a center, center and saddle 
point. Thus we compute the following index sum 

indE + indE* + indQ 3 + indQx + indQ 2 + indQ 4 + indQ 5 + indQ 6 = 2(-l) + 5(+l) + (-1) = 2, 

which is consistent with the Poincare-Hopf Index theorem. 

Depending on the algebraic criteria (involving the vortex strengths), there are exceptional cases where Q5 
merges with Q\ or Q4 coincides with Q2, or possibly both, in which case one or both of Q\ and Q2 can assume 
the form of a degenerate isolated stationary point of index +2. In all cases, the corresponding index sum 
can be shown to be in agreement with the Poincare-Hopf formula. Exceptional cases notwithstanding, once 
again we also have, for any parameter values consistent with the hyperbolic case, a complete characterization 
of three vortex dynamics owing to the unique integral path lifting principle. 

Parabolic case (^ 2 ^= 0) 

Three vortex dynamics represented in the trilinear coordinate phase plane for an example of the parabolic 
case is illustrated in Fig. 3. Observe, as was discussed above, in this case the points E and E* lie on the 
singular curve C in T and its copy C* on T* comprised of stationary points of the system (13). Unfortunately, 
the unique integral path lifting principle does not apply on the double D (C) of C, which is a circle on M 
to which the points E and E* both belong, and with this is associated another problem with the singular 
curve exhibited in Fig. 3; namely, the trajectories of (13) appear to cross this curve, which is not what one 
expects for a reasonable dynamical system in which standard uniqueness theorems obtain. But the principle 
does apply on the complement of D (C) on M, making it possible to obtain an almost complete description 
of three vortex dynamics in the parabolic case. To fill in the gaps, it is necessary only to conduct a more 
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Figure 2: Hyperbolic type with K\ = 2, k 2 = 1, and k 3 = — 4/5. 
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careful analysis of the dynamics of the full system in a neighborhood of the singular circle D (C) as in Ting 
et aZ.|5rj]. but we shall not pursue this further here since it not needed for our analysis in the sequel. 

Note that if we cut the sphere along the singular curve D (C) in the example shown in Fig. 3, we obtain 
two disks B\ and B2 for which we can apply the extension of the Poincare-Hopf formula to manifolds with 
boundaries. The upper disk B\ contains the stationary points Qx, Q2 and Q§ in its interior, and it is easy 
to verify that these fixed points are, respectively, a center, center, and saddle point. Whence, we find that 
the index sum in B\ is 

2(+l) + (-l) = l = X ( J B 1 ). 

For the lower hemisphere B2, we find that Q3 is the only interior fixed point, and it is a center of index +1. 
Accordingly we compute that, as expected, 

l(+l) = l = x(B 2 ). 

Depending on the parameter values of the system for a particular example of the parabolic case, certain 
dynamic properties can vary from the example depicted in Fig. 3. However, as shown in Tavantzis & Ting 
1551 (and can also be inferred from the Poincare-Hopf formula in virtue of the fact that Q3 is the only 
stationary point below the singular curve C), the fixed point Q3 is always a center. 

Before leaving our summary of the trilinear coordinate phase plane behavior of the three vortex problem, 
we wish to emphasize the fact that regardless of the type, the stationary point Q3 is always a center, and 
there are often additional centers. 

III. TYPES OF PERTURBATIONS 

Here we shall provide examples of perturbations three point dynamics that satisfy the hypotheses of the 
main results that we shall derive in succeeding sections. These particular types of perturbations will also 
be used for numerical illustrations of where the theorems apply, and where they begin to break down as 
evidenced by the onset of nonregular, chaotic regimes. The types of perturbations that we choose are three 
vortex dynamics in a half-plane, restricted four vortex dynamics, and three coaxial vortex ring dynamics in 
space. 

A. Three vortex dynamics in a half-plane 

The complex Hamiltonian governing equations for three point vortices of nonzero strengths Tk at respec- 
tive points Zk — Xk+iyk, 1 < < 3, in motion in an ideal fluid in the half-plane H := {z = x+iy 6 C : y > 0} 
are 

3 3 
Kkh = i K j K k(zj - ZkT 1 - i^KjKkizj - z k y x =2id Zk H A , (1 < k < 3) (21) 

where the Hamiltonian function is 

Ha ■= - ^ K i Kk lo S \ z i ~ z k\+ X] K i Kk lo S 1% ~ Zk \ ■ ( 22 ) 
i<i<fc<3 i<j<fc<3 

With the understanding, expressed in Section II, that smoothness requirements actually demand that the 
domain be restricted to points in := {z € H 3 : Zj 7^ z^ \fj 7^ k}, we can express (21) in smooth, real 
Hamiltonian form 

x k = Kl l d Vk H A = {H A ,x k } , Vk = -^ l d Xk H A = {H A ,y k } , (1 < k < 3) (23) 

where we use the same Poisson bracket as in (4). This system, as is easily verified, has the following 
independent constants of motion in involution 

Ha, 3 J ■= Kiyi + n 2 y2 + ^3^3, (24) 

and this appears to be the maximal such set of invariants: Although we are unaware of a proof showing that 
(21), or certain analogs for vortex dynamics on a portion of a sphere such as in Kidambi & Newton |35j . is 
in general not LA-integrable, careful numerical studies such as that of Knio et al. |36| provide compelling 
evidence of the existence of chaos for the three vortex problem in the half-plane. 
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From a perturbation perspective, we can obviously write the Hamiltonian function in the form 

H A =H +H AA , (25) 

where 

Has ■= ^ K 3 K klog\zj - z k \ . (26) 
i<i<fc<3 

Now it is obvious from the nature of these various functions that there exist extensive regions in H# at any 
positive distance away from the boundary (the x-axis) where we have both \H a \ -C |i?o| and \d z H a \ <C \d z Ho\, 
in keeping with our point of view of treating (21) as a (small) perturbation of three vortex dynamics. 

B. Restricted four vortex dynamics 

By the restricted four vortex problem (dynamics), we mean the motion of four point vortices in an ideal 
fluid in the (complex) plane C (— K 2 ), where three of the vortices, located at points Z\, z-i and z 3 , have 
respective nonzero strengths Ti, T2 and T3, while the fourth vortex at the point Z4 has strength satisfying 
I I <C m := min{|ri| , | T2 1 , 1^1} and we have the freedom of making |T4| /m as small as we wish. It follows 
from (4) that the complex Hamiltonian governing equations are 

4 

K k z k = i KjK k (zj ■- z k ) _1 = 2id Zk H B , (1 < k < 4) (27) 

where the Hamiltonian function is 

H B := - ^ KjK k log\zj - z k \ ■ (28) 

\<j<k<i 

Of course this can be put into the following real Hamiltonian form using the same Poisson brackets as in (4): 

Xk = K^dy^B = {H B , x k } , y k = -K^d Xk H B = {H B ,y k } , (1 < k < 4) (29) 

where analogous adjustments in the domain are obviously required to insure smoothness. Just as in the case 
of (4), symmetry considerations show that (26) has the motion invariants 

4 4 

H B , J := ^KkZk, K := J^Kfc \z k \ 2 , (30) 

k=l k=l 

with the following independent integrals in involution 

H B , K, L := (SHJ) 2 + (3J) 2 . (31) 

These three independent integrals in involution represent the maximal number, as it has been proved that 
the system (29) is not integrable in general (see Bagrets & Bagrets [H] and Ziglin [1H|)- We note, however, 
that there are certain special cases of four vortex motion on a plane or sphere that are LA-integrable, as 
shown in Aref & Stremler 0, Borisov et al. |17j . Eckhardt |23| . Sakajo |49| . and Sokolovskiy & Verron |5()| . 

Our focus here is to treat (27) or (29) as a perturbation of the three vortex problem. To be mathematically 
precise, we actually have to consider the system as a perturbation of the three vortex problem embedded in 
C 4 , by which we mean 

3 

K k Zk = i KjKk (zj ■- z k y x — 2id Zk H a , (1 < k < 3) 

l 4 = = 2id Zi H . (32) 

Obviously this embedded system is also LA-integrable, and its dynamics consists of a two-parameter infinity 
of copies of three vortex dynamics - one for each constant value of 2:4. This having been said, there is no 
harm in the slight abuse of notation that we shall use from now on of treating the restricted four vortex 
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problem as a perturbation of the three vortex problem. To highlight the perturbation aspect, we write the 
Hamiltonian function in the form 

Hb = H o + Hb,i, (33) 

where 

3 

H B ,i ■= -cr^Kj log \zj - z 4 \ , (34) 

3=1 

and we have replaced k 4 by a to emphasize the fact that this parameter may be chosen to be as small as 
necessary to insure the desired perturbation behavior. 

C. Three slender coaxial vortex ring dynamics 

Consider three (circular) coaxial vortex rings of respective nonzero strengths T\, T2 and T3, with axis 
of symmetry the y-axis, in motion in an ideal fluid in K 3 . The rings intersect any meridian half-plane 
bounded along the y-axis in unique points (ri,yi), (r '2,2/2) and (^,2/3), respectively, where r is the distance 
from the y-axis, and the motion of the rings is completely determined by these points owing to the axial 
symmetry of the equations of motion. It is convenient to set z := r 2 + iy — s + iy, so that the motion of 
the three intersection points can be considered to be in the half-plane H. The complex Hamiltonian form of 
the dynamical equation of motion of these points - de-singularized in a standard way to eliminate the usual 
infinity in the self-induced velocity of the rings -can be expressed as (cf. Blackmore & Knio HHU, Blackmore 
et al. JI], and Lamb |38j) 



" pit/ z 

z k = V" KjKkrjrk (zj - z k ) / ^7*- 3/2 cos2a ^ a ~~ iK t r k 1 [ lo S ( 8r fe^J 1 ) _ °- 55 ^ 

3 „tt/2 

— 2i KjKkTj / (rj — rk cos 2a) A 

„•_-, .-ji. Jo 



da (35) 



2id Zk H c , (1 < k < 3) 



where S c is a very small positive number representing the common core radii (used in the de-singularization 
procedure) of the vortex rings, 

A ifc := (r fe - rjf + (y k - y,) 2 + Arjr h sin 2 a, (36) 
and the Hamiltonian function is 

3 „tt/2 

H C := -ti^iejrj [log (Sr^" 1 ) - 1.558] - 4 ^ ^K^n / AJ^ /2 cos2ada. (37) 

j=l l<j<k<3 Jq 

This system can, employing the same Poisson bracket used in (4) , be recast in the real Hamiltonian form 
Sk = Kk 1 d yk Hc = {H c , Sk }, yk = -^ 1 d Sk H c = {Hc,y k }- (1 < k < 3) (38) 
It is easy to show that (37) has the two following independent integrals in involution 

3 3 

iir.c-. »vv ».<•':;• ( 39 ) 

3=1 i=i 

However, there are no additional independent invariants in involution, as proved in Bagrets & Bagrets jS], so 
although the dynamics of two coaxial rings is LA-integrable, this is not true in general for the system (35) 
or (38). 

The formulation of (35) as a perturbation of three vortex dynamics is rather more subtle than that of the 
two preceding examples. Formulas obtained by Callegari & Ting |18j indicate that, relative to a coordinate 
system moving along the axis of symmetry with the overall translation velocity of the ring configuration, the 
equations of motion of the rings are closely approximated (modulo a constant factor) by those of three point 
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vortices at the intersection points of the rings in a meridian plane, when the rings are sufficiently close to 
one another compared to their distance from the axis of symmetry. More precisely, Blackmore & Knio (Hj 
showed that in a coordinate system moving with the center of vorticity 

z cv := KlZl+K2Z2 + K3Z3 =J/fiM, (40) 

which is defined for any planar configuration of point vortices as long as ^ 0, the equations of motion 
are Hamiltonian, with a Hamiltonian function of the form 

H c = - (5KW2) [Ho + H x ] , (41) 

where \H\\ / \Hq\ = o (1/ logp) and |9 z i?i| / \d z Ho\ — o (p) as p — > 0, where p is the ratio of the diameter of 
the configuration of meridian plane points of intersection of the rings to the distance of z cv from the axis 
of symmetry. As 9tz cl , is a constant of motion of (35), the above analysis establishes the three coaxial ring 
problem as a perturbation of the three vortex problem for rings. We note that the above formulas require 
that ^ 0, and, to simplify matters, we shall usually assume hereafter that this is the case. One can 
prove the results that we obtain in the sequel without this assumption, but certain rather straightforward 
modifications are required in our methods of proof in some instances. 



IV. QUALITATIVE REGULARITY RESULTS 

To get an idea of the type of results that we shall present in this section, we first state a theorem for 
the examples of the preceding section for the cases in which all of the vortex strengths have the same sign. 
The proof of this theorem is either contained in, or follows directly from, the results in Blackmore & Knio 
Blackmore et al. 11 , Blackmore et al. Blackmore & Champanerkar Khanin [32], and Ting 
& Blackmore 55 . 

Theorem 1. Suppose that in each of the perturbation examples in the preceding section the parameters Kj 
all have the same sign, so that ^ and the center of vorticity is defined. Then the following properties 
hold with respect to a moving coordinate system with origin at the center of vorticity: 

(i) For the three vortex problem in the half-plane (23) governed by the Hamiltonian function Ha there 
is a set of initial configurations of positive (Lebesgue) measure with the diameter of the configuration 
sufficiently small with respect to the distance from the x-axis such that the motion is quasiperiodic on 
invariant KAM tori, interspersed with periodic orbits. 

(ii) The restricted four vortex dynamics governed by Hb according to (29) exhibits, for a set of initial 
conditions of positive measure, quasiperiodicity on invariant KAM tori along with periodic orbits. This 
holds when the diameter of the initial configuration is sufficiently small and a is sufficiently small, or 
the fourth vortex is sufficiently far from the three larger vortices, or for some combination of both the 
conditions on a and the distance of the smaller vortex from the larger vortices. 

(hi) Three coaxial vortex ring dynamics generated by He in accordance with (38) exhibits the following 
behavior: There is a set of positive measure of initial configurations with the distances among them 
sufficiently small compared to their minimum distance to the axis of symmetry that produces quasiperi- 
odic motion on invariant KAM tori and periodic orbits. 



Proof. The proof, as indicated above, follows directly from the cited papers. □ 

In the remainder of this section, we shall find rather general conditions on a perturbation term H\ of a 
Hamiltonian system 

t k = {if, z k } , (42) 

where l</c<3, orl</c<4in the embedded version described above for the restricted four vortex 
problem. Here the Poisson bracket is the same as used in all of the previous equations, the Hamiltonian 
function is of the form 

H = H + H U (43) 
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and the assumptions are general enough to subsume the persistence of regularity results of Theorem 1, as 
well as including the same types of perturbations when the vortex strengths associated to the three vortex 
dynamics generated by Hq are allowed to differ in sign. 

To accommodate all of the perturbations discussed in Section III, we shall consider H to be defined and 
analytic on an open subset of an appropriate complex unitary space, such as for the half-plane and 
coaxial ring problems, and C 4 ^ for the restricted four vortex problem. Moreover, we shall also include the 
possibility of the perturbation H\ depending on a real (freely chosen) parameter a such that 

\H 1 /H \,\d m H 1 \/\d m H Q \^0 (44) 

uniformly as |er| — * on any compact subset of the domain on which these expressions are defined. It is 
convenient to introduce some additional notation. Wc shall refer to a subset of the domain of definition 
of (42) as ample if it is of positive (Lebesgue) measure in the domain or a submanifold of the domain of 
complex dimension at least two. Moreover, we shall say that the system (or H) has the uniform perturbation 
property ( UPP) if there exists a sequence S — {S n : n 6 N} of subsets of the domain of the system with 
nonempty interiors such that S n +i C S n , 

\Hi/H \ < / and \dM / \d z H \ < — j- (45) 

log(n +1) n+ 1 

for every z £ S n for all neN, where N denotes the positive integers. In this case, we call the sequence S 
a uniform perturbation filtration (UPF) for (42). We are now in a position to state the main result of this 
section in a very concise way. 

Theorem 2. Let system (4-2) satisfy the UPP as described in (45), and also (44) when the perturbation 
depends on a parameter as in the case of the restricted four vortex problem. Furthermore, suppose that for 
some UPF, S = {S n : n S N}, the system satisfies the property that there is a companion sequence of compact 
subsets of positive measure JC = {K n : n € N} with K n C S n for all n £ N such that 

z(0) e K n+1 =^ z(t) e S n (46) 

for all t > and for each sufficiently large n, where z{t) is the solution of (42) initially at z(0). Then there 
is an ample set of initial conditions for the system (42), and a set of sufficiently small values of \a\, when 
there is such a free parameter dependence as in the case of the restricted four vortex problem, for which the 
dynamics includes quasiperiodic motion on an ample set of invariant KAM tori interspersed with a countable 
collection of periodic orbits. 

Proof. Our argument relics heavily on what might be called the limit KAM theorem and a useful generaliza- 
tion of the Poincare-Birkhoff theorem, which were successfully employed by Blackmorc and his collaborators 
|H1 1111 IXrtj to prove the existence of, respectively, an ample set of invariant tori and periodic orbits for several 
examples of vortex dynamics problems of the type under consideration. 

First, following the same approach as in those papers, we rewrite the Hamiltonian in terms of the action- 
angle coordinates associated to the unperturbed LA-integrable system as 

W(A,9) = Wo(A)+Wi(A,e), (47) 

where the action and angle vectors are, respectively, A = (Ai, A/.) and = (9i, 0/.), with 2 < k < 4. 
We have used slightly different notation for the various Hamiltonians to underscore the fact that it may be 
necessary, as in the case of coaxial vortex rings, to reduce the number of degrees of freedom by one in order 
to insure the nondegeneracy of Tlo, which naturally alters the forms of the original terms of the Hamiltonian 
(43). This also explains why we have indicated that the number of degrees of freedom k in (47) can assume 
values between two and four. More specifically, k = 2 for the coaxial vortex ring problem, k = 3 for three 
vortex in a half-plane dynamics, and k = 4 for the restricted four vortex problem (in which case we recall 
that the three vortex problem must be embedded in C 4 ). 

For our purposes it is not necessary to know the exact form of H.q (A), which is very difficult to deduce 
in general; we need only take note of the following readily verifiable properties that hold in any compact 
subset of the domain in which (42) is analytic: 
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/ d 2 n (A) dHp(A) \ 

A 2 := det y |^ °£ J * 0, (49) 

and 

W (A)/A 1) A 1 /A 2 ^0 (50) 

as |A| 2 :=A? + ... + A2 ^0 {cf. 0). 

Translating the hypotheses to the transformed Hamiltonian system generated by Ti. (A, 0); namely 

A, = d @j H = d @j Hi, &j = -d Aj H = -d Aj H - d Aj Hi, (51) 

we infer the existence for any e > of a compact set K e of positive measure of initial conditions for (51) 
such that 

\Hx (A(t), 0(f)) /Ho (A(t))| , |3( A ,e)Wi (A(t), 0(t))| / |d A W (A(t))| < e (52) 

for all t > 0. Consequently, owing to the nondegeneracy condition (48), together with (50) in the limiting 
case for |A| — > , it follows from the KAM theorem that there exists an ample set of invariant (real) k- 
dimensional tori for (42). If k = 2, the isoenergetic nondegeneracy condition (49), coupled with (50) when 
|A| — > 0, is enough to insure the existence of periodic orbits. For k > 2, the existence of cycles for (42) can 
be proven, for example, by using the generalization of the Poincare-Birkhoff fixed point theorem employed 
by Blackmore et al. Thus, the proof is complete. □ 

It is a relatively straightforward task to demonstrate that the three types of perturbations described in 
the preceding section satisfy the above hypotheses, thereby demonstrating the existence of ample regular 
dynamics regimes (like those for LA-integrable systems) in each of those cases; to wit, the following result 
is a direct consequence of Theorem 2. 

Corollary 1. The three vortex in a half-plane system (21), which we denote as (A), the restricted four vortex 
problem (27), identified as (B), and the three coaxial vortex ring system (35), denoted as (C), all satisfy the 
hypotheses of Theorem 2, so each of these examples exhibits ample sets of initial configurations leading to 
quasiperiodic flows on invariant KAM tori accompanied by periodic orbits. More specifically, admissible 
UPF's of initial conditions leading to regular motion for each of these examples can be described as follows 
for the elliptic (ET), hyperbolic (HT), and parabolic (PT) types of the unperturbed system: 

(ET) For (A) there are four kinds of UPF's , denoted as S^ A , S A ,S^ A and S A \ with the following 

characterizations: S A ~* is comprised of approximately equilateral configurations ( associated to the point 
E or E* described in Section II) of the three vortices with diameter sufficiently small compared with 
their distance from the x-axis; S A consists of configurations with \Z\ — Zi\ \Z3 — — \Z3 — Zx\ 
(associated with Q3 defined in Section II) of diameter sufficiently small compared with the distance 
of the configuration from the x-axis; S\ , associated with Q2, is characterized by making the obvious 
changes in the definition of S^> an d > corresponding to Q\, is defined by making the evident 
revisions of the description of \ There are also four kinds of UPF's for (B), <Sg , \ Sg ^ and 

Sg\ defined analogously to , S A \S A 2 ^ and S A \ respectively, wherein the ratio of the diameter 
of the configuration of the main three vortices to the distance from the x-axis is replaced by the ratio of 
the diameter to the distance of the configuration from the fourth small vortex; and four kinds for (C), 
Sq \ Sq \ Sq^ and Sq\ also defined analogously to S A *^ , S A ,S A and S A \ respectively, where the 
point vortices are replaced by the intersection points of the rings with any given meridian plane, and 
the distance from the axis of symmetry replaces the distance from the x-axis. 

(HT) There is at least one kind of UPF for (A) corresponding to S A ^ as described above, and as many as 

four others associated with S A ^ and S A \ and possibly two others corresponding to the points Q4 and 
Q5 when they are centers, as discussed in Section II. For (B) and (C) there are analogous UPF's that 
always include those associated to and respectively. 
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(PT) For (A), (B) and (C) there is always one UPF corresponding to 



and Sq\ respectively, 



possibly additional UPF associated to any other possible centers as described in the trilinear phase plane 
characterization of the dynamics covered in Section II. 



Proof. We shall provide detailed arguments only for the three vortex system in the half-plane, either when 
the strengths of all vortices have the same sign (the elliptic type), or they differ in sign (which includes both 
the hyperbolic and parabolic types). The verifications for the restricted four vortex, and three coaxial vortex 
ring problems for the various types of the unperturbed (three vortex) system can be obtained analogously by 
straightforward - but rather lengthy - calculations, and we shall leave the details to the reader, noting that 
the desired results when the unperturbed system is of elliptic type actually follow directly from Theorem 1. 

Even though, as mentioned above, the conclusions we are seeking arc a direct consequence of Theorem 
1 when the unperturbed system is of elliptic type, we shall include a proof here, since it will be helpful in 
pointing the way to establishing sufficient conditions for the hyperbolic and parabolic types. Recalling that 
there is no loss of generality in assuming that K\ > n 2 > n 3 > when the unperturbed system is elliptic, we 
first rewrite the constants of motion (24) of (23) in a form more useful to our purposes; namely 



H A = ci = k\ log (2yi) + k\ log (2y 2 ) + n\ log (2y 3 ) - k y k 2 log 



{xi - x 2 ) 2 + (yi - y 2 f 



(x 1 - x 2 ) 2 + (y 1 +y 2 ) 2 



1/2 



log 



(xi - x 3 ) 2 + (j/i - y 3 ) 2 



(xi -x 3 ) 2 + (yi + y 3 ) 2 
3 J = c 2 = K\yi + n 2 y 2 + n 3 y 3 . 



1/2 



- K 2 K 3 log 



(x 2 - x 3 ) 2 + (y 2 - y 3 ) 2 



(x 2 - x 3 ) 2 + (y 2 + y 3 ) 2 



1/2 



(53) 



and define M (c l7 c 2 ) to be the set of points in satisfying the pair of equations (53). 

It is straightforward to show directly from the form of these defining equations that given any arbitrarily 
large and small positive number, respectively, A and v, there exist A* = A* (A, v) > A and < = v*{\, v) < 
v such that the following property is satisfied if c\ and c 2 are chosen to be sufficiently large positive numbers: 
let the initial positions of the vortices, z\ (0) = xi(0) + iyi(0), z 2 (0) = x 2 (0) + iy 2 (0), z 3 (0) — x 3 (0) +iy 3 (0) be 
approximately in the shape of an equilateral triangle (corresponding to E in Fig.l) of diameter less than f*, 
with yi(0), 2/2(0), 2/3(0) > A*, and define M* (ci, c 2 ) to be the component of M (ci, c 2 ) containing the initial 
point z(0). Then for all configurations (z\,z 2 ,z 3 ) in the component M* (c\,c 2 ), the diameter is less than or 
equal to v, and 2/1, 2/2)2/3 > A. Owing to the connectedness of orbits of (23), it must therefore follow that the 
diameter of the configuration z(t) = (z\(t) , z 2 (t) , z 3 (t)) is less than or equal to v, and yi(t), y 2 (t),y 3 (t) > A 
for all t > 0, where z(t) is the solution of (23) with the specified initial condition. Whence the construction 
of the desired UPF associated to these sets is a simple matter. 

If the unperturbed system is hyperbolic or parabolic, then as indicated in Section II, we may - and do - 
assume without loss of generality that K\ > n 2 > > n 3 . This puts a slightly different complexion on the 
system of equations (53), which must be satisfied by all solutions of (23). However, not so different that 
we cannot use the same kind of argument as for the elliptic type system (given above) modulo a few rather 
obvious modifications. In fact, we can take our cue for the necessary adjustments by recalling our discussion 
of the trilinear phase portraits for the three vortex problem in Section II. It is not difficult to see from a 
close inspection of (53) for the case when k 3 is negative - in which we concentrate on those terms containing 
the negative vortex strength - that we can simply change the initial condition on the vortex configuration 
to conform to the center Q 3 (see Figs 1 and 2). More precisely, we need only change the description of the 
initial configuration in the preceding paragraph so that \z\(0) — z 2 (0)\ <C |^i(0) — 2:3(0) | ~ 1 2:2(0) — 2:3(0) | 
in order to obtain orbits of (23) staying in a component of M(ci,c 2 ) analogous to M* (01,02). Just as in 
the previous paragraph, this leads directly to the desired UPF's for the hyperbolic and parabolic cases, and 
completes the proof for three vortex dynamics in the half-plane. □ 

Before moving on to a study of perturbations of specific periodic orbits of three vortex dynamics, we note 
that the following generalization of Theorem 2 can be easily proved by using essentially the same arguments 
as in its proof given above. 

Theorem 3. Suppose the system (4-2) is a perturbation of a general LA-integrable Hamiltonian system 
generated by Hq, and that H is defined and analytic on an open subset of C m , with m > 1. In addition, 
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assume that (42) satisfies the hypotheses of Theorem 2, and that the usual KAM nondegeneracy condition 
holds for Hq in an admissible UPF. Then there exists an ample set of initial conditions (and small values of a 
parameter, if pertinent) such that the system exhibits quasiperiodic flows on an ample collection of invariant 
tori, together with periodic orbits. 

V. PERSISTENCE OF PERIODIC ORBITS 

In the preceding section we proved theorems demonstrating that various types of, generally non-integrable, 
Hamiltonian perturbations of certain LA-integrable systems - including that governing the motion of three 
point vortices in an ideal fluid in the complex plane - have ample dynamical regimes exhibiting the regularity 
properties that characterize integrable systems; namely quasiperiodic motion on invariant tori and periodic 
orbits. 

Such qualitative results as those already obtained, which are essentially existence theorems, beg the 
question of how closely such behaviors of the perturbed system approximate those of the unperturbed LA- 
integrable system associated to Hq? This can be viewed as a more quantitative classical perturbation theory 
related query about the qualitative entities whose existence has been proven by more modern methods in 
symplectic dynamics. In this section we provide a partial answer to this question as it relates to periodic 
orbits, which is embodied in the following result. 

Theorem 4. Let the Hamiltonian system defined by (42) and (43) be a perturbation of three vortex dynamics 
generated by Hq with respect to a coordinate system moving with the center of vorticity. Suppose that the 
closed curve C represents a periodic orbit of the unperturbed three vortex system with Hamiltonian function 
Hq. Define the (compact) tubular neighborhood X e (C) for each positive e as 

1 e (C) :={z: A(z,C) <e}, 

where A denotes the usual (unitary) distance function. Also define 

K{C) = m^{\H x /H \ , \d*H x \ I \d z H \ : z e T £ (C)} . 

Then if for a given e small enough to insure that % e (C) is a smooth ( = C°° ) submanifold of the domain 
in which the system (42) is smooth, the quantity A £ (C) is sufficiently small, (42) has a periodic orbit C in 

T e (C). 

Proof. First select a point zo € C, and let £q and £ be the energy hypersurface, respectively, of the 
unperturbed system generated by Hq and the perturbed system generated by H, which contain the point 
zo. Both £q and £ have a common odd real dimension, which we denote as 2m — 1. Let A r denote the 
Ricmannian metric in E induced by the metric A, and define 

B s :={ze£: A r (z,z ) <<$} 

for small positive values of 5. This is clearly a (2m — l)-ball in £ for all 5 < Sq sufficiently small. 

We let (fit represent the flow generated by (42). The following properties follow directly from standard 
results on differential equations (see e.g. when A e (C) is chosen to be sufficiently small: There exist 

a transversal S through Zo for the system (42) in £ and < S^ < Si < So such that (i) each of the sets 
(3g := Bg n X is a 2(m — l)-ball in £ whenever < 6 < 6%, (ii) the orbits of (42) initially on j3g 2 all 
pass through the interior of (3g t in finite (positive time) generating a Poincare map P : f3g 2 — > f}$ 1 with 
P (f3$ 2 ) C interior(/3,5 2 ), (iii) the radial geodesic curves in E emanating from ip t (zo) all transversely intersect 
the boundary of tpt(/3g 2 ) in unique points for all < t < t m , where t rn is the maximum first return time to S 
among all those for the flow of (3s 1 , and (iv) the flow of [3s 1 generated by (42) through the first return time 
to E is contained in T e (C). 

Our proof is obviously complete if P has a fixed point, so assume on the contrary that this is not the 
case. As [3s 1 and [3$ 2 are both (real) odd-dimensional balls, it follows from properties (i)-(iv) that we can 
modify (42) in the flow of f3s 1 \Ps 2 in order to obtain an extension P of P, such that P maps (3g 1 into itself 
and has no fixed points in /3g 1 \f3$ 2 . Thus, this smooth self-mapping P of /3g 1 has no fixed points at all, 
which contradicts the Brouwer fixed point theorem. Accordingly we conclude that P must itself have a fixed 
point, which corresponds to a cycle C of (42) in T e (C), so the proof is complete. □ 
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Considering the generality of the above result, and the relative simplicity of the proof, it seems as though 
this should have certainly been discovered. However, the authors were unable to find such a result in the 
literature, although it appears that it could be derived rather directly from certain results that prove the 
existence of periodic solutions for Hamiltonian systems employing index theory (c/. Blackmore & Wang jlU) . 
Gole and Josellis PO] )• It is interesting to take note of the very special case where d z Hi in the above 
theorem takes the form d z H\ = (f>(z,)d z Ho, where is a smooth, real valued function whose magnitude can 
be made arbitrarily small in the tubular neighborhood T e (C). Then as one can readily see by making the 
obvious transformation of the time parameter, the cycle of the perturbed system is not just an approximation 
of C, it is identical with it, and the local flow of the perturbed system is identical with that of the unperturbed 
system (modulo a change of parametrization). 

Theorem 4 can be applied directly to the perturbations of three vortex dynamics that we have been 
considering to identify conditions under which the perturbation has a periodic orbit close to one for the 
three vortex problem. We leave the very straightforward proof of the next result to the reader. 

Corollary 2. The hypothesis and conclusions of Theorem 4 regarding the existence of a periodic orbit of the 
perturbed system close to a periodic orbit C : z — z(t), < t < p, of the unperturbed (three vortex) system 
hold for (A) three vortices in a half-plane, (B) restricted four vortex dynamics, and (C) three slender coaxial 
ring dynamics obtain under the following conditions: 

(a) Each of the coordinates of z(t) on C is sufficiently distant from the x-axis in the complex plane. 

(b) For the coordinates z(t) — (zi(t), z 2 (t), z 3 (t), Zi(t)) G C, the distance in the complex plane C from 
24(0) to the set {z G C : z — zi(t),Z2(t) or 23(2;), < t < p} is sufficiently large, or the magnitude 
of the strength of the fourth vortex is sufficiently small, or a suitable combination of both of these 
conditions is enforced. 

(c) The coordinates of z(t) on C representing the points of intersection of the rings with a meridian plane 
are sufficiently far from the axis of symmetry (represented by the y-axis in this plane). 

We note that the conditions given in Corollary 2 by no means exhaust all possible situations where the 
perturbed system has a periodic orbit close to a periodic orbit of the unperturbed system. For example, it is 
easy to see that in the restricted four vortex problem there are many such cases where the fourth small vortex 
is neither particularly small nor very distant from the three large vortices: Simply consider any configuration 
of the large vortices that generates a periodic solution of the three vortex problem, and place the fourth 
vortex, of any strength, at the center of vorticity of the larger vortices. Then the fourth vortex remains 
fixed, and the motion of the three larger vortices is unaffected by its presence. 

VI. SIMULATIONS OF DYNAMICS 

In this section we provide numerical examples that illustrate the persistence of regular motion for different 
types of perturbations, as well as breakdown of regularity as large perturbations are considered and the 
hypotheses of our theoretical results are accordingly violated. As mentioned earlier, we consider half-plane 
dynamics, restricted four vortex perturbations, and slender coaxial vortex ring dynamics , which we refer to 
below as type A, B, and C, respectively. 

Our simulations cover just a small sample of possible cases for the various perturbations, and they are 
presented in two basic types of graphical forms: As trajectories in the plane or half-plane, or as Poincare 
sections, which are especially well suited to illuminating transitions from regular to chaotic motion. In 
particular, for the unperturbed system we present plots of the trajectories of the three vortices in the plane, 
juxtaposed with the corresponding trilinear phase plane stationary point. If any of the perturbation types 
satisfy the hypotheses of Corollaries 1 and 2, we expect that plots of the trajectories of their vortex elements 
(with respect to a coordinate system moving with the configuration) would show a small variation of the 
plots for the unperturbed system. Also included for type A perturbations are Poincare maps, along with 
the corresponding Poincare map for the unperturbed system, showing that there is a transition to chaotic 
motion as the array starts closer to the boundary of the half-plane compared to the diameter of the initial 
configuration of vortices. For type C perturbations, we present Poincare maps exhibiting strong regularity 
when the slender coaxial rings are very close to one another compared to the distance of the configuration 
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from the axis of symmetry, as expected in view of Corollary 1. For the restricted four vortex problem, a 
pair of Poincare maps shows how the motion tends to be more chaotic as the initial position and strength 
of the fourth small vortex, respectively, starts closer to the initial group of three larger vortices and grows 
in comparison to the strengths of this group. This behavior is entirely consistent with Corollaries 1 and 2. 

For type A and C perturbations, we consider the setup illustrated in Fig. As shown in the figure, the 
geometry is specified in terms of two parameters, b and a, while the strength of the perturbation is reflected 
by the height h for type A, or the mean radius p for type C. Introducing the parameter c = v 'a 2 + 6 2 , an 
analysis was conducted of the nine cases summarized in Tabled 



a k 



h or p 



Figure 4: Initial configuration for type I and II perturbations. 



Case 


b 


c 




K2 


^3 




1 


0.1 


1 


2 


1 


0.5 


3.5 


2 


0.1 


1 


2 


1 


-2/3 





3 


0.1 


1 


2 


1 


-0.8 


-0.4 


4 


2/7 


1 


2 


1 


0.5 


3.5 


5 


2/7 


1 


2 


1 


-2/3 





G 


2/7 


1 


2 


1 


-0.8 


-0.4 


7 


0.49 


1 


2 


1 


0.5 


3.5 


8 


0.49 


1 


2 


1 


-2/3 





9 


0.49 


1 


2 


1 


-0.8 


-0.4 



Table 1: Summary of inputs. The calculated value of is also reported. 



Summary of Simulation Results 

• Our main focus is on examples where the unperturbed system is of the elliptic type, i.e. cases 1, 4, 7. 
Initial configuration near Q3, in between two separatrices (region 3B shown in Fig. 1), and near the 
center E. We expect more sensitivity to perturbations for of any of the three types for case 4 than 
1 and 7, since case 4 represents a configuration closer to the separatrix shown in Fig. 1, and this is 
borne out by our simulations. 

• The motion of the three vortices in the plane for cases 1, 4, and 7 in the unperturbed system is shown 
in Figure [5] Notice that even in this case, the trajectories are considerably more complicated for case 
4 than cases 1 and 7. Again this is due to proximity to the separatrix. 
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For a type A perturbation with h = 10, we observe little impact on trajectories of the 3 vortices. This 
is consistent our theoretical results. 

For cases 1 and 7, we also consider lower values of h - as low as h = 2.5. For case 4 corresponding, to a 
point close to the separatrix for the unperturbed system, the dynamics become much more complicated 
at a more rapid pace as h decreases. The Poincare maps (x%,yi for yi — 0) shown in Fig. [5] suggest 
that when h = 2.5 the dynamics is starting to exhibit chaotic regimes. 

Also considered in our simulations are type C perturbations. We found that it is necessary to start 
with a large mean radius for the three slender coaxial vortex rings in order to observe theoretical 
predictions. With p — 250, the motion remains regular. This is illustrated in the Poincare maps shown 
in Fig.d 

In the final set of simulations, we considered examples of restricted four vortex dynamics The configu- 
ration of the three larger (principal) vortices, which we designate as vortices 1-3, is defined in terms of 
equal-strength vortices, k% = 1 for i = 1, 2, 3, initially located on the vertices of an equilateral triangle 
of side 1. Specifically, the initial conditions are given by: zl = (—0.5,0), z2 = (0.5,0), and 23 = (0,h) 
where h = y/3/2. For the unperturbed case, the trajectories with respect to the centroid are identical 
circles centered at the origin. We consider perturbation due to the introduction of a weak vortex having 
&4 = 0.1. Five different initial conditions for the fourth vortex are considered, namely Z4 = (0.1,0), 
Z4 = (0.1, h/3), Z4 = (0.1,2/i/3), Z4 — (0.1, h), and Z4 = (0.1, Ah/3). We refer to these as cases i - v, 
respectively. 

Consistent with predictions, the trajectories of vortices 1-3 are only weakly affected by the introduction 
of the fourth (weaker or smaller) vortex (not shown) . More precisely, we find that the common circle of 
motion of the three principal vortices becomes slightly thickened (forming a thin annulus) in the plane 
- with the thickness of the annulus and variation of motion within it reflecting the magnitude (as well 
as the particular form) of the perturbation caused by the smaller vortex. One expects the annulus to 
thicken when the smaller vortex starts closer to the configuration of vortices 1-3, its strength increases, 
or when the weaker vortex starts in a position that has a more subtle connection with the configuration 
of the whole set of four vortices. These subtleties will be discussed at some length in the following 
points. 

The most dramatic effect on the dynamics for the restricted four vortex problem is naturally manifested 
in the motion of the fourth (smaller or weaker) vortex. The trajectory of the fourth vortex can be 
either regular or complex and even chaotic, depending on its initial location and strength. Moreover, 
the complexity of its trajectory is quite sensitive to changes in its initial location and strength. This 
is illustrated in the Poincare maps of (xi , j/4 at y\ = 0) . 

The Poincare maps for the motion of the fourth vortex shown in Fig. [S] indicate via the characteristic 
splattering that the dynamics for cases i and ii exhibit typical chaotic behavior, or at least a transition 
from regular to chaotic motion, whereas cases iii and iv are essentially regular. This may appear to 
be somewhat surprising in light of the fact that the distance from vortex 4 to the set of vortices 1-3 in 
cases i and ii is not very different from the distances in cases iii and iv. 

A plausible explanation for the somewhat counter intuitive behavior described above is as follows: If we 
view vortex 4 as perturbing the motion of vortices 1-3 (see Fig. 1), we would not, as indicated above, 
expect this to have much effect on the dynamics of vortices 1-3. In particular, the configuration of 
vortices 1-3 places it initially at the center E, so the inherent stability should preserve periodic motion 
for this configuration that varies only slightly from an equilateral array. Dually, we can consider the 
triple of vortices 1, 2 and 4, as being perturbed by the strong third vortex. For cases i and ii, the 
initial configuration of {1,2,4} is rather close to a separatrix shown in Fig. El This experiences a strong 
perturbation from vortex 3, which is substantial enough to push the configuration across the separatrix 
and also possibly break this curve. Similarly, cases iii-v are such that the initial configuration of vortices 
i, ii and iv are substantially further removed from such a separatrix, so then it is not surprising that 
the dynamics for cases i and ii is far less regular than that for cases iii and iv. 
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• Some of these subtle points are further illustrated and contrasted in Fig. EI which shows the dynamics 
of the fourth vortex, again using Poincare maps. The dynamics of the fourth vortex of finite strength, 
portrayed in the left column of figures, is contrasted with the same cases on the right, where the 
strength of vortex 4 is taken as zero. In other words, we are treating vortex 4 as a passive particle 
propelled by the motion of the three stronger vortices. Notice that in all cases where the strength of 
the vortex is zero, the Poincare maps show that the motion of this vortex is regular. The juxtaposition 
in Fig. [5] underscores the point that the effect of the fourth vortex is dependent both on its initial 
placement and strength. 
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-0.5 -0.3 -0.1 0.1 0.3 0.5 




-0.5 -0.3 -0.1 0.1 0.3 0.5 




-0.5 -0.3 -0.1 0.1 0.3 0.5 



Figure 5: Trajectory of the three vortices with respect to the centroid (left) and in the a — f3 plane (right): 
case 1 (top), case 4 (middle) and case 7 (bottom). The trajectory of vortex 1 is depicted in black, while 
those of vortices 2 and 3 are shown in red and green, respectively. 
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Figure 6: Poincare maps for case 4. Left: unperturbed system; middle: type A perturbation with h = 10; 
right: type A perturbation with h = 2.5. 

0.4 | . , , , . , , 1 

0.3 - 




-0.3 - 



-0.4 -0.2 0.2 0.4 

Figure 7: Poincare maps for case 1 (black), 4 (red), and 7 (green). Type C perturbation with p = 250. 




-2-10 1 2-2-10 1 2 



Figure 8: Poincare maps for the restricted 4- vortex configuration. Left: maps for case i (black), iii (red) and 
v (green). Right: maps for case ii (black) and iv (green). 

VII. CONCLUDING REMARKS 

Following a brief summary of the dynamics of three point vortices moving in an ideal fluid in the plane 
in both a Hamiltonian and trilinear coordinate context, which highlighted the properties most pertinent for 
our investigation - such as the existence of a center regardless whether the unperturbed system is of elliptic, 
hyperbolic or parabolic type corresponding, respectively, to K1K2 + + K2K3 positive, negative or zero 
(assuming that Ki + K2 + K3 7^ 0) - we described the three kinds of perturbations that comprised the focus 
of this paper. These were three vortices in a half-plane, a restricted four vortex problem, and the motion of 
three coaxial vortex rings. 

We then formulated and proved some new results on the existence of regular regimes - characterized by 
behavior associated to integrable Hamiltonian systems - for non-integrable perturbations of LA-integrable 
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Figure 9: Poincare maps for the restricted 4-vortex configuration, for perturbed (left) and unperturbed 
(right) systems. Shown are results for cases i, ii, hi, and iv, arranged from top. 



three vortex dynamics. The regular dynamics in point includes the existence of an ample set of quasiperiodic 
flows on invariant KAM tori along with periodic solutions. These results, which were proved employing rather 
modern variants and generalizations of the KAM and Poincare-Birkhoff fixed point theorems, extend some 
related recent results for point and coaxial ring vortices to cases where the signs of the vortex strengths 
differ, and are such that their sufficient conditions can be readily checked. But, they are really just existence 
theorems that provide little information about the location and finer geometric details of the predicted 
invariant tori and periodic orbits. 

The lack of quantitative specificity in our first set of perturbation theorems was subsequently addressed 
for the case of periodic orbits, with some novel results concerning the existence of periodic solutions for 
perturbed systems that closely approximate cycles of the unperturbed three vortex system. These results 
- possessing the look and feel of classical perturbation theorems - were proved using Poincare sections in 
concert with the Brouwer fixed point theorem, are stated in terms of easily verifiable criteria guaranteeing 
the existence of periodic orbits of the perturbed system close to those of the three vortex system. 

Several numerical simulations were performed to illustrate the perturbations theorems obtained in a 
variety of contexts. Also included were examples to indicate how the desired regularity of the dynamics 
breaks down as the assumptions of our theorems are stretched to the limit, which is accomplished by showing 
how there is a (transitional) preponderance of chaotic regimes as these limits are exceeded. 

Our study of perturbations of the three vortex problem has, we hope, shed some new light on these 
particular systems as well as more general vortex phenomena. It certainly has - as is often the case in such 
studies - created more interesting questions than it has answered, several of which we intend to investigate 
in the near future. For example, our work here suggests that most of our results can be extended to systems 
comprised of a larger number of vortex elements if they arc clustered in certain ways. As an illustration 
of this, suppose that one has a pair of vortices of strengths of the same sign that are initially quite close 
together. Then this pair constitutes a binary system that should, under the right conditions, remain close 
together for all time, and behave as if it were a single vortex of strength equal to the sum of the strengths of 
each of its components if any other point vortices have initial distances from this pair that are substantially 
larger than the distance between the pair. A configuration of four vortices including such a binary system 
and two other vortices initially quite distant from the pair, might well be expected to behave very much like 
a slightly perturbed three vortex system in the plane (with the binary system acting like a single vortex), 
thus generating a profusion of regimes exhibiting regular (integrable type) dynamics. Interesting clustering 
problems of this kind are something that we plan to turn our attention to in our future research. 

Another rather manifest question raised by our investigation - that we intend to presently turn our at- 
tention to - concerns extending the quantitative perturbation results obtained for periodic orbits to invariant 
tori. It appears that it may be possible, using some type of generalization of the Poincare-Birkhoff fixed 
point theorem, to show that in certain cases the perturbed system has an invariant torus that is close to an 
invariant torus for the unperturbed system. 
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